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Abstract 



Quenched lattice QCD simulations and quenched chiral perturbation theory are used 
together for this study of strangeness in the nucleon. Dependences of the matrix elements 
on strange quark mass, valence quark mass and momentum transfer are discussed in both the 
lattice and chiral frameworks. The combined results of this study are in good agreement with 
existing experimental data and predictions are made for upcoming experiments. Possible 
future refinements of the theoretical method are suggested. 
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I. INTRODUCTION 

The effects of virtual strange quarks on the properties of a single nucleon represent basic 
information about QCD and the strong interaction. Hence, there is presently a great deal 
of enthusiasm for studies of the nucleon's strangeness electric and magnetic form factors. 
Recent experiments have produced two measurements]!], § and ongoing efforts are expected 
to provide more results soon[[|. 

First principles calculation from QCD requires the use of lattice field theory techniques, 
and a number of explorations have been carried out by various authors @, ||, ||, |7], . The pres- 
ence of the disconnected strange quark loop and the smallness of the resulting strangeness 
form factors cause lattice simulations to be expensive and the extraction of meaningful 
results to be difficult [0]. 

Chiral perturbation theory (ChPT) can play a valuable complementary role alongside 
lattice QCD. ChPT is QCD's low-energy effective theory written in terms of the physi- 
cal hadrons rather than quarks and gluons, and it contains low energy constants (LEC's) 
whose numerical values should be determined from lattice QCD or directly from experi- 
ment. Quenched SU(3) ChPT|], [T(J] corresponds to quenched QCD with three active quark 
flavours — up, down and strange — and it produces analytic expressions for the strangeness 
form factors that explicitly display their dependences on the strange quark mass, valence 
quark mass and momentum transfer. It is clearly advantageous to relegate as much of the 
calculation as possible to ChPT so that valuable computer time can be spent on the physics 
that ChPT cannot predict. In other words, one need only extract the required LEC's from 
lattice QCD simulations, and then the strangeness form factors can be studied directly in 
quenched SU(3) ChPT. 

On the other hand, the strangeness form factors can in principle be measured in lattice 
QCD simulations with minimal recourse to ChPT: the strange quark mass and the momen- 
tum transfer can be fixed to their physical values in a lattice simulation and then ChPT 
is only needed for extrapolation of the valence quark mass. This extrapolation can be per- 
formed with quenched SU(2) ChPT rather than SU(3), thereby providing the benefit of a 
more rapid convergence for the chiral expansion since it no longer requires expansion in 



powers of the strange quark mass ||1 1|| . 



In the present work, we report the results of high-statistics lattice QCD simulations 
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for the strangeness electric and magnetic form factors together with the strangeness scalar 
density. A number of different analysis methods are employed and found to give consistent 
results. Two strange quark masses, three valence quark masses and five momentum transfer 
values are studied. We also present the analytic quenched SU(3) ChPT formulae for the 
three strangeness matrix elements of interest and apply them to our lattice QCD data. 
The alternative of using quenched SU(2) ChPT is briefly discussed as well. Finally, we 
compare our results to the existing experimental measurements, make predictions relevant 
to upcoming experiments, and suggest directions for future theoretical work. 

Our main conclusions are that the raw lattice results for the strangeness electric and mag- 
netic form factors (before any use of ChPT) are very small, that ChPT-based extrapolation 
to the physical up and down quark mass region does not substantially change this, and that 
the lattice QCD predictions are therefore consistent with existing experimental results. 



II. NUMERICAL SIMULATIONS 



The gauge field configurations used in this study were generated from the Wilson gauge 
action at (3 = 6 on 20 3 x 32 lattices, corresponding to a lattice spacing of 

a = 0.1011(7) fm (1) 

as obtained by the authors of Ref. |I2| from a physical string tension of \/K = 427 MeV. 
Actually the lattice spacing is not uniquely determined in the quenched approximation, and 
the authors of Ref. || used the physical nucleon mass to arrive at a = 0.115(6) fm. Our 
full ensemble of 2000 configurations was produced from various independently thermalized 
Markov chains. Within each chain, either 2000 or 5000 triple-step heatbath updates (i.e. 
applied to three SU(2) subgroups) were executed between saved configurations. 

The Wilson fermion action was used to obtain three valence quark propagators per con- 
figuration, having n v = 0.152, 0.153 and 0.154. These correspond to pion masses of 

am T = 0.4772(9)^(13], 0.4237(8)p| and 0.364(1)13 (2) 

respectively. The valence quarks in our simulations have Dirichlet time boundaries; the 
source is four timesteps away from the boundary. On our 20 3 x 32 lattices, the five smallest 
momentum squared values are 

a 2 q 2 = n(7r/10) 2 , n = 0, 1,2,3,4. (3) 
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Tabulated in Table [I] are the energies of a nucleon having degenerate quarks and each of 
these momentum squared values. 

Strangeness matrix elements are calculated using standard methods. This involves a 
three-point function in which a strange-quark loop is correlated with the nucleon propagator. 
It is prohibitively expensive to compute the strange quark loop exactly at every lattice site, 
so we employ a stochastic estimator with real Z 2 noise[|T^J. To reduce the variance, the 
first four terms in the Ki expansion («j denoting the loop quark's hopping parameter) of the 
quark matrix were subtracted for the strangeness electric and magnetic form factors, and 
the first five terms were subtracted for the scalar density |i6||. This stochastic estimation 
method is unbiased with any number of noises, and the statistical uncertainties associated 
with this noisy estimator decrease as one increases the number of noises and/or the number 
of gauge field configurations. 

For K\ = 0.152 we have computed a 60- noise estimate for each of our 2000 configurations, 
and for k\ = 0.154 we have computed a 200-noise estimate for 250 configurations. The vector 
meson masses for these Ki values are 912(8) MeV and 1066(4) MeV respectively (see Table 
VI of Ref. ||T2"| ) which surround = 1019 MeV so that our data will allow interpolation to 
a strange quark loop. From the lattice simulations, three ratios are constructed, 

Rx ^^- G (2) (t> d)G(*)(tf,g) ' ( ) 

where Rs, Rm and Re correspond to the scalar, magnetic and electric cases respectively, 

t is the sink timestep and t' the current insertion timestep. The two-point and three-point 

correlators are shown diagramatically in Fig. |l|. 

Strangeness matrix elements are extracted from the ratios of Eq. (|J) . Denoting the matrix 

elements by M with an obvious subscript, these are related to form factors by 

M {s ,m,e}M = l^\^^,Gf) . (5) 

In the magnetic case, i, j and k run over spatial directions and the corresponding indices 
on M M are suppressed for notational simplicity. 

There are various ways in which the matrix element can be extracted from the ratio. For 
example, one can sum the contributions for the strange quark inserted at different times t' . 
One wayflHfl to do this is 

t 

Rxit, t', q) -> constant + tM x (t, q). (6) 

t'=i 
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A disadvantage of this kind of method is that the matrix element does not emerge directly. 
A fit to the time dependence, which in practice may be linear only over a limited range, is 
required to determine Mx- For this reason we prefer a differential method 

t+i 

[Rx(t, t', q) -R x {t- 1, q)} - M x (t, q) (7) 

t'=l 

which gives Mx directly. For completeness we also consider the relation 

^ fixed 

Rx{t, t', q) -> constant + tM x (t, q), with t fixed > t (8) 

t'=i 

used in Ref. ||. 

Finally one has to relate the lattice matrix element to the continuum one. The physi- 
cal scalar density requires wavefunction renormalization and we use the tadpole-improved 
factor, pi 

(N\Ss\N) = (i-P)g%\ (9) 



4n r 



with k c = 0.157096(28)ig 3 []r3|]. The conserved vector current was used for Gj$ and G%\ 
and its normalization is such that no wavefunction renormalization factor is required. 

Fig. ^| shows our lattice data for the scalar density versus timestep, with k v = 0.154 and 
K\ = 0.152, analyzed using Eq. fl7|). In this case there is a very clear signal and, for each value 
of the momentum transfer, the plateau begins about ten timesteps from the source, although 
uncertainties grow with q 2 . Figs. |3| and f| show the magnetic and electric data from Eq. (^) 
with the same k v , Ki values. In contrast to the scalar density, there is no apparent nonzero 
signal. However, using the scalar density results, which suggest that the plateau region 
begins about ten timesteps from the source, as a guide, one concludes that the form factors 
are consistent with zero within uncertainties less than 0.1 for all q 2 values studied. We have 
verified that Eqs. (0) and @ produce compatible results for all three matrix elements. 

The results of fitting each of our lattice measurements to Eq. (|7|) over four consecutive 
timesteps, beginning ten timesteps from the source in every case, are tabulated in Table |Tl 
with statistical uncertainties obtained from a bootstrap analysis employing 3000 bootstrap 
ensembles. If the uncertainties simply scaled with the square root of the number of config- 
urations then the ratio of uncertainties between k v = 0.154 and 0.152 should be near 2.8, 
but the increased number of noises per configuration for k v = 0.154 could reduce this ratio. 
According to Table [0], on by G^ shows a noticeable dependence on the number of noises. 
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These results for G M can be compared to the findings of Ref . || , since those authors also 
work with the Wilson action with the same (3 and k values, although their lattice volume is 
smaller. From 100 configurations with 300 complex Z 2 noises analyzed using the method of 
Eq. (pD only, those authors interpreted their results to imply a nonzero value for G^ . Our 
studies (see Ref. [0] for a specific discussion) suggest that a clearer picture is attained with 
a larger sample of gauge configurations. According to Table [HJ even the small statistical 
uncertainties of the present work do not permit a definitive nonzero determination of G^ . 
The same is true for Gg . 



III. CHIRAL EXTRAPOLATIONS 



Consider quenched SU(3) ChPT with explicit fields for the pseudoscalar meson octet 
(M), spin-1/2 baryon octet (B), spin-3/2 baryon decuplet (T) and external electromagnetic 
and scalar fields. The ChPT Lagrangian is 



*- _ ^-M "r '-'MB "r ^MB ' *-MB ^MB *-MT '"MBT ■ 



■ j 



(10) 



where a superscript "(n)" denotes an nth order contribution from the expansion in the 
smaller scales — momentum transfer, meson masses and the T-B mass splitting A — rela- 
tive to the larger scales A x « and baryon masses. The leading loop diagrams for our 
three strangeness form factors begin at third order and are displayed in Fig. |5|. Each dia- 
gram receives contributions from various quark flows which have been calculated using the 



approach of Labrenz and SharpefllO]- Besides these loop contributions, there are also contact 
terms in the Lagrangian which contribute low energy constants (LEC's) to the strangeness 
matrix elements. Here are the explicit formulae: 



(N\ss\N) = d/i + C£(A)/iA- 
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G^V) = C 3 + CX(A)A + 
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where A > and 



A(x) 



y/ z — A 2 arccos (;^=J 



for A < y/z, 
for A > ^/i, 



with 



(14) 



(15) 



z = m^- — x(l — x)g . 

Our interest is in spacelike q 2 , so z is positive definite throughout the range < x < 1. The 
g 2 of each lattice data point is obtained from 

2 



q 2 = (E n - E ) 



n 



\10aJ 



(16) 



where 77 = 0, 1, 2, 3 or 4 and the E n are taken from Table |. 

Cb contains the familiar axial couplings (D and F) and Ct contains the octet-decuplet 
coupling C (defined, for example, in Ref. |T0|): 

5 



C B 



^D 2 — 2DF + 3F 2 , 
C 2 . 



(17) 
(18) 



The parameters C±, C2, ■ ■ ■ C% are LEC's, some of which depend on the dimensional regu- 
larization scale A such that the full matrix elements are independent of A. 7 is the ChPT 
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parameter for the quenched r/jnj and m Ss is the mass of a doubly-strange pseudoscalar 
meson. The normalization convention corresponds to F n ~ 93 MeV and \i is the ChPT 
parameter defined by 

m K = /i(m + m s ) (19) 

with m = m u = m^. In order to verify various aspects of these ChPT expressions for 
(N \ss\ N) (q 2 ), G^(q 2 ) and Gg\q 2 ), comparisons were made to the collection of papers in 
Ref. |9|. 

Notice that the three strangeness matrix elements contain a total of six parameters - 
fiCi + /iAC£(/i), C 3 + AQ(/i), C 5 + Cg(/i)/mjv, Cb, Ct and 7 2 — and the dependences 
on each of these parameters are linear. Because D, F, C and 7 are real parameters it 
follows that Cb, Ct and 7 2 must be positive definite, and Eq. (|TT| ) therefore requires that 
(N \ss\ N) (q 2 ) decreases as m, m s or — q 2 is increased. This is consistent with the lattice 
QCD data of Table y. 

It should be noted that the range of q used in our lattice simulations extends far beyond 
the range of applicability of ChPT, and there is therefore no reason to expect that the form 
of ChPT will look anything like the lattice data for these larger momentum values. As 
would be hoped, use of only the lattice data at smaller momentum values leads to a good 
ChPT fit. As it happens, the ChPT expressions fit all three matrix elements surprisingly 
well over the entire momentum range studied. Although this is surely accidental, it means 
that the ChPT expressions can be used as a convenient method of smoothly interpolating 
the momentum dependences of these matrix elements. 

To determine numerical values for the six parameters appearing in the ChPT expressions, 
we perform a least squares fit to the data of Table [TT| In particular, we'll fit the 39 data 
points having k\ = 0.152 (data for G^(0) are omitted since gauge invariance requires a zero 
result) and verify that predictions for Ki = 0.154 are consistent with our lattice simulations. 
We will also perform an independent fit using only 12 of the 39 data points: those having 
a 2 (f = or a 2 (f = (71-/IO) 2 . These smallest momenta are the ones most appropriate to 
ChPT and, as will be demonstrated, the final predictions for strangeness matrix elements 
are rather insensitive to whether or not the higher momentum data are used as input for 
the ChPT fit. The statistical uncertainties of the fit parameters are determined from a 
bootstrap analysis. 

In addition to the statistical error there is a systematic uncertainty due to the choice of 
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chiral model. The dynamics of the ChPT expressions reside in the loop diagrams, and they 
contain the quenched 7/ parameter r y 2 as well as the non-7/ parameters Cb and Ct- It is 
possible to obtain a good fit to the k\ = 0.152 data in the extreme limit of no rf (y 2 = 0) or 
in the opposite extreme of "maximal r/" where Cb = Ct = 0. (In the maximal 77' case, we 
also choose C3 + C£(A)A = since it is clear from Eq. (|T^) that this parameter would simply 
be an additive constant for G^.)(q 2 ) and would be consistent with zero when fitted to our 
lattice QCD data.) These separate possibilities indicate that our lattice data are not precise 
enough to determine the fraction of 77' physics in the strangeness form factors. One might 
expect the physical values for these parameters to lie somewhere between the two extremes, 
and we will use this range to define a theoretical error bar. The results of our fits to the 
Ki = 0.152 data, and the resulting predictions for k\ = 0.154, are recorded in Table |TTT| . The 
fits are consistent with the direct lattice QCD simulations of Table |I[ The corresponding 
ChPT parameter values are listed in Table [IV], along with the parameter values obtained 
from fits to the data having the two smallest momenta: a 2 (f = and a 2 (f = (7r/10) 2 . In the 
unquenched theory r y 2 does not appear and standard phenomenology leads to Cb ~ 0.9 and 
1.4 < Ct ~ 2. Not surprisingly, the quenched parameter values in Table [IV] are different 
but are still 0(1). 

For physical meson masses, 

777 m 2 1 



, (20) 
m s 2m K — 25 

which leads to 



m Ss = m K \ : , ; , — = 1.39m*:. (21) 
V 1 + m/m s 

At — q 2 = 0, Gg vanishes identically. Figs. |6| and |7| show the other two strangeness matrix 
elements as functions of the kaon mass. Fixing m* to its physical value leads to the mo- 
mentum dependent strangeness matrix elements of Figs. § and which are our final results. 
Comparison to experiment, along with disclaimers about such a comparison, are contained 
in Section fV[ 

To conclude this section we return to the suggestion from Ref. |ll] of using SU(2) ChPT 
instead of SU(3). This is an appealing idea because SU(2) ChPT typically converges more 
rapidly. In effect, the kaon loop diagrams of Fig. [5] get replaced by SU(2) LEC's. Although 
SU(3) ChPT uses a common set of parameters (Cb, Ct and •y 2 ) for the kaon loop effects 
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in all three strangeness matrix elements, SU(2) ChPT has separate LEC's for each ma- 
trix element. Since the raw lattice QCD data of Table [TT] only reveal a nonzero signal for 
the strangeness scalar density, it is difficult to discuss SU(2) ChPT extrapolations of the 
strangeness electromagnetic form factors in any detail. Perhaps future lattice QCD data for 
these form factors will be precise enough to benefit from SU(2) ChPT. 



IV. DISCUSSION 

The results of this work (Figs. ^| and |9|) compare favourably to the available experimental 
data: 

M _ f 0.14 ±0.29 ±0.31, Ref. M, 
G ( m{qI) = U (22) 

[ 0.05 ± 0.06, This work, 

M _ M _ f 0.025 ±0.020 ±0.014, Ref. 01, 

G^(g 2 2 )±0.39GS(g 2 2 ) = ^ (23) 

[ 0.07 ±0.05, This work, 

where — q\ = 0.1 GeV 2 and —q\ = 0.477 GeV 2 . Here, the uncertainties (incorporating 
both statistical and theoretical modeling errors) in our results have been estimated by the 
requirement that all curves from Figs. [8] and |9|, representing fits to all momenta, fits to 
only small momenta, "maximal //" fits and "no 77"' fits are within one standard deviation 
of the quoted central value. The lack of a fundamental scalar probe makes the strangeness 
scalar density harder to extract from experiment, but Figs. |8| and |9| can be compared to 
other quenched lattice QCD simulations. The renormalization group invariant quantity 
representing the fractional strange quark contribution to the nucleon mass is: 



m s (N\ss\N) (0) 



m N 



0.302(48) at p = 5.7 Ref. [20], 

0.195(9) at = 6.0 Ref. 0, (24) 

0.21(11) at @ = 6.0 This work. 



If the curves of Fig. || are not included in the predictions for these strangeness matrix 
elements and if the statistical errors of Fig. |8] are ignored relative to the theoretical errors 
(reflecting the difference between "maximal rf" and "no 7/" fits), then one arrives at the 
earlier results reported in Ref. §: G^(qf) = 0.03 ± 0.03, G%\ql) + 0.39G$(gf) = 0.027 ± 
0.016 and (m s /m N ) (N \ss\ N) (0) = 0.15(2). 
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There are a number of ways that future theoretical studies could improve upon the 
results obtained in this work. From the outset we have restricted ourselves to the quenched 
approximation, and this introduces a systematic error that is perhaps 10-20%p2]|. It is also 
not obvious that higher orders in the ChPT expansion are small for the case at hand, i.e. 
SU(3) ChPT for baryons with quark masses in the strange region. It would be interesting 
to see the results of partially quenched simulations and lighter valence quarks for these 
strangeness matrix elements. Refinements of the disconnected loop techniques could also 
be advantageous, such as perturbative subtraction beyond 0(k 4 ,k 5 ) and heatbath noise 



methodsf23|1. Finally, we recall that the so-called strangeness electric and magnetic form 
factors would not be exactly zero even in a world without any strange quark, due to isospin 
violation [p4| , [25[| . Based on Ref. ||25|| , the isospin violation effects are not so different in 



magnitude from the tiny strange quark effects discussed in the present work. 

Although there are certainly further steps that can be taken toward a more detailed 
understanding of these strangeness matrix elements, the present study has established that 
G^\q 2 ) and G$(q 2 ) are small over the range of momenta and quark masses used in these 
lattice QCD simulations, and that they remain small when extrapolated with quenched 
SU(3) ChPT in combination with lattice QCD data for (N \ss\ N) (q 2 ). 
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TABLE I: Dimensionless energies of the nucleon, aE n , with momentum squared a 2 q 2 = n(-7r/10) 2 . 
All fits begin 16 timesteps from the source, 2000 configurations are used, and statistical uncertain- 
ties are from a bootstrap analysis with 3000 bootstrap ensembles. 



n k = 0.152 k = 0.153 k = 0.154 

0.869(2) 0.799(2) 0.728(3) 

1 0.927(3) 0.862(3) 0.795(4) 

2 0.986(4) 0.924(5) 0.865(7) 

3 1.034(7) 0.977(10) 0.922(15) 

4 1.070(12) 1.013(18) 0.945(30) 
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TABLE II: Fits to the matrix elements of Eq. (^) beginning 10 timesteps from the source. The 
momentum squared is a 2 q 2 = n(7r/10) 2 . Statistical uncertainties are from a bootstrap analysis 
with 3000 bootstrap ensembles. 



Kl = 0.152 



ki = 0.154 



G 



0.152 2.6 

1 1.7 

2 1.2 

3 1.1 

4 0.7 
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TABLE III: Predictions of quenched SU(3) ChPT, after a least squares fit to 39 lattice data. 
The estimated uncertainties include the range between the two extreme cases of maximizing or 
minimizing the quenched rj' contribution in ChPT loop diagrams relative to non-7/ physics, as 
discussed in the text, as well as the statistical uncertainties from a bootstrap analysis. 
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TABLE IV: The parameter values obtained for the two extreme fits to our quenched lattice QCD 
data at k\ = 0.152, as discussed in the text, using (i) lattice data from all available momenta 
and (ii) lattice data with aPq 2 = and aPq® = (7r/10) 2 only. Statistical uncertainties are from a 
bootstrap analysis with 3000 bootstrap ensembles. 
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FIG. 1: Two-point and three-point correlators that appear in Rx of Eq. (||). Each solid line 
represents a quark propagator, and the shaded box denotes a current insertion. 




FIG. 2: Lattice data for the strangeness scalar density as obtained from Eq. (fj]) with k v = 0.154 
and ki = 0.152. Panels (a) through (e) correspond to momenta a 2 q 2 = n(7r/10) 2 with n = 
through 4 respectively. Uncertainties are calculated from 3000 bootstrap ensembles. 




FIG. 3: Lattice data for the strangeness magnetic form factor as obtained from Eq. (0) with 
k v = 0.154 and K[ = 0.152. Panels (a) through (d) correspond to momenta a 2 q 2 = n(ir/10) 2 with 
n = 1 through 4 respectively. Uncertainties are calculated from 3000 bootstrap ensembles. 
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FIG. 4: Lattice data for the strangeness electric form factor as obtained from Eq. (0) with k v = 
0.154 and ki = 0.152. Panels (a) through (e) correspond to momenta a 2 q 2 = n(-7r/10) 2 with n = 
through 4 respectively. Uncertainties are calculated from 3000 bootstrap ensembles. 
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FIG. 5: Leading loop diagrams for the strangeness matrix elements from quenched chiral pertur- 
bation theory. Dashed, solid and double lines denote octet mesons, octet baryons and decuplet 
baryons respectively. A shaded box denotes a current insertion. 
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FIG. 6: Strangeness matrix elements at — q 2 = as functions of mx- The two solid curves 
represent the extreme cases of maximizing or minimizing the quenched rf contributions in ChPT 
loop diagrams relative to non-?/ physics. ChPT parameters are obtained from a fit to 39 lattice 
QCD data points as discussed in the text, and the thickness of a hatched band denotes statistical 
uncertainties from 3000 bootstrap ensembles. 
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FIG. 7: Strangeness matrix elements at — q 2 = as functions of mx- The two solid curves 
represent the extreme cases of maximizing or minimizing the quenched rf contributions in ChPT 
loop diagrams relative to non-?/ physics. ChPT parameters are obtained from a fit to 12 small- 
momentum lattice QCD data points (a 2 (f = and a?if = (7r/10) 2 ) as discussed in the text, and 
the thickness of a hatched band denotes statistical uncertainties from 3000 bootstrap ensembles. 
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FIG. 8: Strangeness matrix elements as functions of — q 2 . The two solid curves represent the 
extreme cases of maximizing or minimizing the quenched rf contributions in ChPT loop diagrams 
relative to non-7/ physics. ChPT parameters are obtained from a fit to 39 lattice QCD data points 
as discussed in the text, and the thickness of a hatched band denotes statistical uncertainties from 
3000 bootstrap ensembles. 
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FIG. 9: Strangeness matrix elements as functions of — q 2 . The two solid curves represent the 
extreme cases of maximizing or minimizing the quenched rf contributions in ChPT loop diagrams 
relative to non-r/ physics. ChPT parameters are obtained from a fit to 12 small-momentum lattice 
QCD data points (aPq 2 = and a 2 ^ 2 = (-zr/lO) 2 ) as discussed in the text, and the thickness of a 
hatched band denotes statistical uncertainties from 3000 bootstrap ensembles. 



